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As  a  preliminary  we  examine  a  two -ray  condition  as  shown  in  Fi,';ure  4  . 
One  wavefront  is  assumed  to  arrive  al< ng  the  x-axis,  another  at  an  angle  0. 
The  wavefronts  are  sinusoidal  in  time  and  the  sum  at  points  along  x  is 

s(t,x)  =  Aj^cos(ut  -  kx)  +  A2Cos^a)t  -  ({>  -  kx  cos9)  (1) 


FIGURE  4.  RAY  ARRIVALS 


Tlie  case  treated  here  Is  analogous  to  nne  encountered  in  HI  systems  with 
sinu-oidal  inter forenct*  to  a  desired  carrier. 


QPR  No.  ?4 


Aj^  and  A2  represent  the  magnitude  of  the  two  wavefronts,  ^  is  the  phase  dif¬ 
ference  between  them  at  x  ■=  0,  and  k  =  ItifX  where  X  is  the  wavelength.  A 
sensor  at  point  x  responding  equally  and  linearly  to  both  wavefronts  would 
see  s(t,x).  Writing  (1)  in  envelope-angle  form. 


where 


s(t,x)  =  B(x)cos(mt  +  l'(x)) 


B(x) 


+  A„ 


2Aj^A2  cos  [(f  -  kx(l  -  cos6)] 


i|)(x)  =  tan 


-1 


Aj^sin  kx  +  A2sin(iji  +  kxcosS) 
Aj^cos  kx  +  A2Cos(i{i  +  kxcosS) 


(2) 

(3) 

(4) 


The  phase  obtained  using  (4)  will  be  inodulo-2T.  It  is  useful  to  deal  with 
the  phase  derivative  dijj/dx  if  the  modulo-2Tr  ambiguity  is  to  be  avoided.  It 
can  be  shown  that 


“  k  -  |(1  -  cos6) 


where 


1  + 


1  -  a 


1  +  a^  +  2a  cos[({)  -  kx(l  -  cos0)  ] 


a  “  Aj^/A2 


(5) 


(6) 


When  ^  is  large,  meaning  that  the  important  part  of  the  received  w’ave  is 
along  the  x-axis 


=  k  -  cos[4i  -  kx(l  -  cos9)],  (7) 

It  fluctuates  sinusoidally  around  k  with  the  fluctuation  amplitude  decreasing 
to  zero  as  a^  goes  to  infinity.  For  ^  small 

^  =  k  cos6  -  ka(l  -  cosG)  cos[(J>  -  kx(l  -  cos0)],  (8) 

again  a  sinusoidal  fluctuation  which  decreases  to  zero  as  ii  approaches  zero. 
Fol  intermediate  values  of  a^  the  fluctviation  of  d-Ji/dx  is  as  shown  in  Figure  5, 
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FIGURE  5 ,  PHASE  DERIVATIVE  VARIATION  WITH 
DISTANCE  IN  THE  TWO-RAY  CASE. 

The  variation  is  periodic  with  period 


1  -  cose 

and  has  peak  excursions  above  and  below  k  given  by 


(9) 


6  =  k(l  -  cos0)/(a  -  1)  (10) 

a  =  -k(l  -  cos9)/(a  +  1)  (11) 


Figure  5  is  drawn  assuming  a  >  1.  The  fluctuation  is  around  k,  which  turns 
out  to  be  the  average  of  dij>/dx.  Note  that  at  a  =  1^  B  is  positive  and  high 
In  magnitude,  and  the  fluctuation  is  highly  impulsive.  For  a  <  1,  di[)/dx 
fluctuates  around  the  value  k(l-cos0)  rather  than  around  k  and  for  a  =  1  6  ir 

negative  and  high  in  magnitude,  and  the  fluctuation  is  again  impulsive  but 
negative  going.  When  ^  is  close  to  unity  the  phase  as  a  function  of  posi¬ 
tion  (which  is  the  integral  of  di()/dx)  is  as  shown  in  Figure  6. 

In  underwater  applications  the  angle  G  typically  found  in  long  range 
paths  is  less  than  20°.  Assuming  it  to  be  10°  the  period  is 


X 


A 

1  -  0.98 


SOX 


(12) 


> 
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It  Is  Interesting  to  note  that  In  analyses  found  in  the  literature  of  the  cor 
relation  distance  of  underwater  acoustic  waves,  numerical  estimates  around 
SOX  are  typically  obtained  (see  for  Instance  II]). 


[1]  P.  Vr.  Smith,  Jr.,  "Spatial  Coherence  in  Multipath  oi  Multlmodel  Channels,” 
Journal  Acoustic  Society  of  Amerlra,  Vol.  60,  No.  2,  Aug.  1975,  pp.  305-310 


We  now  turn  to  the  more  general  case  of  n  wavefronts  arriving  at 
angles  6^,  j  *■  l,2,....n.  The  received  sum  at  a  position  x  is 


s(t,x)  ••  Z  A.cos(i>)t  -  -  kx  COS0.) 

J=1  J  ^  J 


(13) 


The  j  =  l,2,...n  are  random  phase  angles  of  each  of  the  wavefronts  on 
arrival  at  the  point  x  =  0.  It  is  convenient  to  write  this  in  the  form 


s(t,x)  *  B(x)  cos[a)t +i|)(x)  ] 


jwt 


where 


“  Re  z(x)e 


"  j(4>.  +  kx  COS0  ) 

!(x)  =  =  I  A  c  J  J 

j=l  J 


(14) 


(15) 


iKx)  is  the  phase  angle  we  will  study  ar.d  as  was  done  before  we  find  the 
phase  derivative 


dip 


Im(-~  -% 


dx  "  'z(x)  dx^ 


From  (15)  we  have 


dz 


~  *=  j  S  A.  k  cosO.  e  ^ 

dx  ■'j.i  J  J 


j(i}i^  +  kx  cos0j) 


so  that  (16)  becomes 


dip 

dx 


lA,  k  COS0.  e 


j(<|)j  +  kx  cosOj) 


Im  j 


1 


J 


EA.  e 

J  ^ 


j  (<}ij  +  kx  cosOj) 


E  E  A  A,  k  cosO.  c 
i  j  ^  J _ ’ 


j[<f.  ~  ‘♦'j  kx(cosOj  -  cor.0^)l 


j  (iji .  +  kx  cosO  . )  I  ? 

E  A  e  ^ 

J 


E  E  k  cosO^  co;;[(j.^.  --  +  kx(cor.O^.  -  cosOj)! 

E  E  A. A  coi;[i{>  -  .  +  kx(eosO,  -  cii-.O.)! 

J  '  J  J  J  ' 


(16) 


(17) 


(18) 
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e 

As  a  rule  0^^  «  ^  ^nd  cose^  =  1 - Where  this  approximation  is  permis¬ 

sible  ve  can  write 


where 


«=  k  -  ^  ^ 

dx  2  11  a^jCx) 

aij(x)  =  A^Aj  cos[((i^  ”  kx(cos0^  -  cosO^)] 

To  retrieve  i{»  we  must  integrate  dij;/dx;  l.e.  we  form 


(19) 


dil'Cx  ) 

-  dx^  =  tp(x)  -  ^J)(0)  (20) 

o  ’'l 


The  integration  will  generate  the  phase  difference  between  the  phase  at  x 
and  the  phase  at  the  origin  of  integration.  From  (19)  we  see  that  one  term 
on  integration  will  te  kx,  the  linear  phase  variation  associated  with  the 
normal  phase  vs.  position  function  of  a  plane  w’ave  along  the  direction  of 
travel  of  tho  wave.  In  beam  forming  with  an  array  of  sensors  along  x  one 
will  subtract  the  phase  progression  kx  if  the  axis  of  the  beam  is  to  be 
colinear  with  the  x  axis.  In  this  case  the  remaining  phase  difference  be¬ 
tween  a  point  X  and  the  origin  is 


-k 

2 


r.  E  0, 


Li 

r.  f  a 

o  i  j 


aij(x^) 


dx. 


(21) 


2 

If  the  approximation  cosOj^^  =  1  -  0^.  /2  is  not  used,  the  renviining  phase 
after  correcting  for  kx  is  given  by  subtracting  k  from  (18)  and  integrat¬ 
ing  over  X. 

Numerical  evaluations  of  the  remaining  phase  difference  have  been  made 
for  a  number  of  cases. One  particular  case  is  shown  in  Figure  7  determined 
assuming  21  equal  amplitude  rays  arriving  at  2°  intervals  from  0  =  -20*’  to 


* 

riograminlng  of  this  copqnitatlon  and  the  one  described  later  giving  array 
pattern,  war  clone  by  I)o  Juan  Ho. 
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FIGURE  7.  PHASE  DIFFERENCE  vs  POSITION 


+20“  with  respect  to  the  horizontal  (see  Figure  4),  each  with  a  phase  angle 
^  randomly  selected  in  the  Interval  (0,27r).  Because  the  average  wavelength 
of  the  various  rays  as  seen  along  x  is  less  than  X,  the  wavelength  along  the 
direction  of  travel  of  the  ray,  there  is  a  linearly-tending  phase  accumula¬ 
tion  with  distance  as  seen  in  Figure  7-  On  top  of  this  accumulation  there 
is  a  random  variation.  The  fluctuation  around  a  straight  line  approximation 
to  the  phase  difference  ranges  around  +3  radians.  Thus  even  if  the  phase 
were  corrected  to  account  for  the  slope  of  the  straight  line  approximation, 
a  +  3  radian  random  error  would  still  be  encountered. 

Figure  7  was  obtained  with  one  randomly  selected  set  of  ray  arrival 
phase  angles.  Additional  examples  will  be  ultimately  computed  for  different 
sets  of  arrival  phase  angles  to  provide  data  suitable  for  obtaining  statisti¬ 
cal  averages.  Other  cases,  including  different  intervals  of  arrival  angle, 
different  ray  amplitudes,  and  different  numbers  of  arriving  arrays  will  also 
be  treated. 

Having  a  sample  function  of  phase  vs.  position, a  logical  next  step  is 
to  determine  the  gain  and  pattern  of  the  random  planar  floating  array  when 
it  is  focused  in  some  azimuthal  direction  using  conventional  beaaformlng, and 
when  the  source  signal  is  propagating  toward  the  array  through  the  multipath 
medium.  As  a  first  step  a  program  was  developed  for  selecting  element  posi¬ 
tions  over  a  circular  area  assuming  a  uniform  distribution  of  clement  posi¬ 
tions. 

If  the  array  is  assumed  confined  to  a  circle  of  radius  P  with  uniform 
distribution  over  the  circle,  the  density  function  in  the  joint  random  vari¬ 
ables  X,  Y,  is 

Pjj  y(x,y)  =  — 2  »  X  +  y  ^  P 

’  Tf  P 

<=  0  ,  elsewhere. 

Transforming  to  polar  coordinates,  (R,5),  we  have 


ifp 

=»  0 


0  r  ^  p 
0  <  £  2  n 

elsewhere 
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The  marginal  densities  in  R  and  0  arc 


Pj^(r)  ^  0  1  1  P 

P 

®=  0  ,  elsewhere 

P  .  0  <  4,  27r 

*•0  ,  elsewhere 


The  random  variables  R  and  4>  are  independent  and  Independent  choices  of  these 
variables  arc  made.  Sample  values  of  <t>  are  obtained  by  a  conventional  cota- 
puter  program  which  selects  sample  values  uniformly  distributed  in  (0,1)  and 
multiplies  these  by  2ir.  Sample  values  of  R  are  obtained  by  picking  a  number 
Z  uniformly  distributed  in  (0,1)  and  forming 


for  then 

°  fr  “  “I  ’  0  <  r  5  P  . 

P 


Finally,  the  pairs  (r,4)  so  obtained  are  converted  back  into  rectangular  co¬ 
ordinates  by 


X  =  r  cos4 
y  -  r  sin^i 


Using  element  positions  so  determined  the  array  pattern  was  next  found. 
The  geometry  of  the  problem  is  sho'.-.Ti  In  Figure  8. 

Assuming  N  elements  distributed  over  the  circle,  copiiased  to  form  a  beam 
along  the  y  axis,  the  array  pattern  is  given  by 

,  N  J[k(x  coni'  +  y  slnO  -  y  )  -  a(d  )] 

A(»)  -  I  I  i  BW  )e  "  ■>  »  B  I 

n=l 


where  x  ,  y  ,  and  d  are  defined  in  Figure.  8,  and  a(d  )  is  a  phase  vs  po- 
nun  n 

sltlon  function  of  the  form  obtained  earlier  and  shown  In  Figure  ?•  B(d  )  is 

n 

the  arrplltud>'  of  the  acoustic  field  at  tlie.  nth  clomonl  .  'llils  quantity  can  be 
obtained  u.sing  the  earlier  analysis,  but  for  our  pviiposcs  now  we  will  assume 
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it  constant  and  set  it  equal  to  unity  for  all  n.  Amplitude  fluctuations  as 
a  rule,  cause  minor  effects  compared  to  phase  fluctuations.  The  phase  sample 
function  of  Figure  7,  called  now  a(x)  is  used  alone  below  to  assess  the  ef¬ 
fect  of  the  multipath  medium.  The  variable  x  in  Figure  7  is  replaced  by  d  , 

n 

with 

d  =  p  -  r  cos  (A  -  i}>) 
n  n  a 

.  2  ,  2.1/2  .  -1  ^n  .. 

=  p  -  (x  -1-  y  )  cos  (tan - t,  ) 

n  n  X 

Q 

Computer  calculations  of  A(<}i),  as  described  above,  were  carried  out  for  two 

cases:  (1)  a(d  )  =0  and  (2)  a(d  )  as  given  by  Figure  7,  and  the  results 

n  n 

are  showTi  in  Figures  9  and  10.  Case  1  is  that  cf  propagation  through  a  trans¬ 
parent  (non-multipath)  medium  while  case  2  is  for  the  particular  multipath 
case  resulting  in  the  phase  function  discussed  above.  Note  that,  the  gains  along 
the  rcaln  beam  in  the  two  cases  are  in  the  ratio  of  about  A. 4  dB  -  a  substan¬ 
tial  factor;  the  sidelobe  structure  is  different  in  detail  but  not  in  general 
characteristics.  These  results,  it  must  be  recognized,  are  based  on  one  set 
of  random  arrival  phases  and  on  on'’  set  of  random  element  positions;  whether 
they  are  ropresentatlvu  remains  to  be  detorinini-d.  Averaging  over  many  sets 
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FIGURE  9 . 


ARRAY  PATTER!,'  OF  Ri\!.T)OM  PLA!;AR  ARRAY  IN  TRANSPARENT 
MEDIUM.  ENT)FIR£  BE^V-I  FORMED  AT  ARIiaJTH  A1,'GLE,  <p  =  90*. 
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FIGURE  10.  ARRAY  PATTFJ'.N  OF  R/O.TOM  R!.A;:AR  ARRAY  ASSU!-1ING  laM.TIPATH 
PROPAGATION.  IRRR-'IRJ!  REAM  FOrOIEU  AT  AZlM'uTil  /AGLE  (}i  »  ! 
'  FOR  TPj'A’SPARENT  MEDIIEI. 

of  arrival  phases  and  element  positions,  as  well  as  carrying,  out  additional 
computations  with  otfier  system  parameters,  remain  to  be  done. 

Fred  Haber 
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